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Using Monte Carlo simulations and the two e g -orbital model for manganites, the stability of 
the CE and A-type antiferromagnetic insulating states is analyzed when quenched disorder in the 
superexchange Jaf between the ti % localized spins and in the on-site energies is introduced. At van- 
ishing or small values of the electron-(Jahn-Teller)phonon coupling, the previously hinted "fragility" 
of these insulating states is studied in detail, focusing on their charge transport properties. This 
fragility is here found to induce a rapid transition from the insulator to a (poor) metallic state upon 
the introduction of disorder. A possible qualitative explanation is presented based on the close 
proximity in energy of ferromagnetic metallic phases, and also on percolative ideas valid at large 
disorder strength. The scenario is compared with previously discussed insulator-to-metal transitions 
in other contexts. It is argued that the effect unveiled here has unique properties that may define 
a new class of giant effects in complex oxides. This particularly severe effect of disorder must be 
present in other materials as well, in cases involving phases that arise as a compromise between very 
different tendencies, as it occurs with striped states in the cuprates. 
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I. INTRODUCTION 

The understanding of transition metal oxides is 
one of the most important open challenges in Con- 
densed Matter Physics. Exotic and complex phe- 
nomena emerge from the correlated nature of carri- 
ers in these materials, such as high critical tempera- 
ture superconductivity and colossal magnetoresistancc 
(CMR). 1 ! 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 In recent years, 
much progress has been made in this context after the 
realization that several competing states, with quite dif- 
ferent properties, are very close in energy in these com- 
pounds. As a consequence, small variations in carrier 
density, temperature, pressure, magnetic fields, and other 
variables often lead to giant responses and an intrinsic 
electronic softness £ Moreover, the effect of quenched dis- 
order appears to be important to understand the CMR 
effect, according to recent theoretical and experimental 
investigations briefly reviewed below. This combination 
of strong correlation and disorder effects, and the simul- 
taneous relevance of several degrees of freedom such as 
spin, charge, orbital, and lattice, lead to a plethora of 
unusual properties, with tremendous potential for func- 
tionalities. 

In this paper, using theoretical models and Monte 
Carlo simulations, the unexpected giant responses of 
some insulating states to the introduction of quenched 
disorder are discussed. Recent investigations by Aliaga 
et alJ& already briefly described a curious "fragility" 
of the insulating charge-ordered and orbital-ordered CE 
state of half-doped manganites. In those previous inves- 
tigations the term "CE glass" was coined, and it was 
observed that the spin correlations at large distances 
characteristic of the CE state rapidly decrease in mag- 



nitude by increasing the strength of the quenched dis- 
order in the links, namely by adding a random compo- 
nent to the superexchange Jaf between the t? s spins de- 
grees of freedom. In independent studies, Motome et 
alH also observed that disorder suppresses charge order 
more strongly than ferromagnetic order using a one-band 
model for Mn oxides (see also Sen et alm^). These con- 
clusions appear to be potentially important since experi- 
mental investigations by Akahoshi et alm^ and Nakajima 
et a/pBi have observed that for a manganite material in 
the clean limit (namely, with minimal quenched disor- 
der) magnetoresistance effects are absent. However, by 
introducing quenched disorder (via a random distribution 
of trivalent ions in the crystal) then a clear CMR effect 
is observed. Thus, experimental evidence suggests that 
quenched disorder is vital to understand the manganites. 
This is in agreement with the early theoretical investiga- 
tions of Burgy et al. using simple "toy" models with two 
states in competition, supplemented by random resistor 
network calculations £221 In those studies, quenched dis- 
order was needed to obtain a CMR effect at low fields. 

In the investigations presented in this manuscript, we 
have studied in considerable detail and using larger lat- 
tices the fragility effect first noticed in Refs. Il^ll9ll20l 
The focus of the investigations is on the half-doped 
regime of manganites, where several previous efforts have 
established the phase diagrams with good accuracy. In 
particular, it is known that a CE phase is clearly sta- 
bilized at intermediates values of Jaf- Of much im- 
portance for the discussion below is that this phase ex- 
ists independently of the value of the electron- (Jahn- 
Teller)phonon coupling A. Such a curious observation 
was understood independently by at least two groups (e. 
g., see Refs. I25ll26|) . Hotta et al~ provided an explana- 
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tion where the well-known zigzag chains of the CE ar- 
rangement arise from an optimization process between 
two competing tendencies: the double exchange mech- 
anism that favors ferromagnetism and the coupling be- 
tween the localized spins Jaf that favors antiferromag- 
netism. This competition leads to an exotic insulator, 
the CE state, in a relatively narrow region of parameter 
space, interpolating between the more conventional ferro- 
magnetic (FM) and antifcrromagnctic (AF) regimes. Fer- 
romagnetic one-dimensional zigzag chains spontaneously 
emerge, with antiferromagnetic coupling between these 
chains. Moreover, the insulating nature of this state 
arises from the unexpected band insulator properties of 
the zigzag chains at n=0.5. This combination of frustra- 
tional and geometrical factors leads to a stable insulating 
CE state. However, with hindsight clearly its nature can 
be anticipated to be quite fragile since its stability is in- 
duced by competition of very different tendencies. 

In our investigations, we have studied the CE state 
in two (2D) and three (3D) dimensions, as well as an- 
other insulating state, the A-type AF state that is stable 
only in 3D. Small amounts of quenched disorder were in- 
troduced and we have analyzed in detail the influence 
of this disorder over the transport properties of the re- 
sulting states, using canonical techniques^ To our sur- 
prise, a rapid transition from an insulator to a poor metal 
was observed. This effect appears robust: it is simi- 
lar in both 2D and 3D, and both for the CE and A- 
typc insulators, to the extent that the electron-phonon 
coupling is small. In other words, it is the "topologi- 
cal" CE insulator unveiled in previous investigations^ 
that shows the strongest fragility toward quenched dis- 
order. This CE state is stabilized neither by a very 
strong nearest-neighbors Coulomb repulsion nor by a 
very strong electron-phonon coupling, but by more subtle 
effects related to frustration and geometry, thus causing 
its fragility. Also, these states appear to be highly sus- 
ceptible to the introduction of magnetic fields, as also 
discussed in the text. 

How exotic is the insulator-metal transition found here 
when analyzed in a more general framework? Anderson 
localization transitions can be induced via the introduc- 
tion of disorder in a metal. In this case, quenched disor- 
der transforms the metal into an insulator, which is the 
opposite of the transition we aim to discuss in this paper. 
Recent investigations on the effects of correlated disorder 
have shown that the metallic state can survive even in 
one dimensional systems^ but still the quenched disor- 
der always favors the insulating tendencies, which does 
not happen in our case. In general, the appearance of a 
metal upon introducing disorder in an insulator is not a 
common phenomenon if carried out at constant electronic 
density. For instance, it is known that insulator to metal 
transitions appear in conventional semiconductors. But 
they occur via doping and the concomitant creation of 
an impurity band. The transmission of charge in disor- 
dered materials sometimes occurs via tunneling between 
large conducting regions^ Doped polyacetylene has a 



similar behavior, although the metallic phase is gener- 
ated via a qualitatively different process^ A similar phe- 
nomenon occurs in other transition metal oxides, such as 
the high temperature superconductors, where chemical 
doping introduces mobile carriers in the system, render- 
ing the compound metallic or even superconducting^ In 
all these materials, carrier doping away from the elec- 
tronic density corresponding to the insulator is needed 
to achieve the metallic state. This is not the case of the 
study reported in this paper that is carried out at a con- 
stant electronic density n. 

Working at fixed n, an insulator to metal transition can 
occur also if uniaxial stress increases the overlaps between 
otherwise localized electronic states. Our analysis is not 
of this variety, since the hopping amplitudes (the model 
parameters which are the most susceptible to pressure) 
are left untouched in our investigations. 

However, there are simple ways in which a disorder- 
induced insulator to metal transition can be envisioned 
and it can occur by the mere broadening of states by dis- 
order. As the disorder strength increases, the gap will 
tend to close. Examples of this behavior were discussed 
by Dobrosavljevic and collaborators^iiS where dynam- 
ical mean-field theory calculations unveiled a metallic 
regime in between Mott insulator and Anderson local- 
ized regions, at values of the disorder strength compara- 
ble to the Hubbard repulsion. However, we believe that 
this simple picture cannot be the full explanation of the 
numerical results reported here since the metallic tenden- 
cies appear even for small values of the disorder strength. 
In the study of Ref. |^ small amounts of disorder do not 
alter substantially the properties of the Mott insulator, 
but the reciprocal occurs: small values of the Hubbard 
U can rapidly transform an insulator to a metal at inter- 
mediate values of the disorder strength. Moreover, the 
subtle physics of zigzag chains is difficult to capture with 
mean- field approximations. An explanation of the CE 
fragility as given in Ref. |2(] for the one-orbital model re- 
sults of Ref. is also based on the simple broadening 
of peaks idea. However, note that the CE state of half- 
doped manganites at vanishing electron-phonon coupling 
seems more fragile than other insulators, and we believe 
that a more sophisticated explanation of its insulator to 
metal transition will be needed. 

Finally, to further motivate these investigations it is 
important to remark that there is evidence suggesting 
that the CE state of the manganites may actually be of 
the exotic "topological" variety discovered by Hotta et 
almH In fact, recent experimental effortsSiM have chal- 
lenged the common folklore that the CE state contains 
ions with sharply-defined charge 3+ and 4+, since no 
evidence of strong charge order was observed. A possi- 
ble explanation of this puzzle is provided in Fig. 10 of 
Ref. [l8| where the charge disproportionation between the 
even and odd sites of the CE state is shown as a function 
of A. It was there observed that the regime of small A is 
the only one compatible with the new experiments. This 
is also in agreement with the well known fact that a fer- 
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romagnctic metallic state must be very close in energy to 
the CE insulating states, since experimentally the latter 
is easily destabilized by a magnetic field into the former. 
At large A, the studies in Ref. showed that all states 
(FM, CE, and AF) are charge ordered, in contradiction 
with experiments. Thus, it is conjectured here that the 
fragility of the CE state may be of key relevance to un- 
derstanding the manganitcs. Recent experimental results 
arc compatible with this conclusion, namely quenched 
disorder affects the CE phase more aggressively than the 
FM stated 

If we accept that the close proximity of states with 
different properties (such as a metal or superconductor 
and an insulator) is the basic cause of the effect, then it 
follows that the giant effect of quenched disorder should 
be present in many other systems as well. For instance, 
the famous stripesS^l in cuprates have recently been 
found to appear at intermediate narrow coupling regimes 
of phenomenological models for the superconductor- 
antifcrromagnetism competitionj22i£i namely they are 
also a compromise between very different states. As a 
consequence, we here conjecture that Cu-oxide stripes 
should also be highly susceptible to the explicit intro- 
duction of disorder. This conclusion can be extended to 
any intermediate narrow phase of transition metal oxides 
that originates from competing tendencies. 

The organization of this manuscript is the following. 
In Section II, the model is described. Since this model 
has been extensively used in previous investigations, the 
discussion is brief. Section III contains the main results, 
namely the influence of link disorder on the transport 
properties of the system. The analysis is carried out for 
two insulating states (the CE and A-type AF phases) 
and in both two and three dimensions. Evidence for the 
metalization of the system upon the introduction of link 
quenched disorder is presented. Section IV describes the 
results for transport now including magnetic fields. Sec- 
tion V deals with on-site disorder, Section VI provides 
evidence that our results are not dramatically altered by 
increasing lattice sizes, and finally conclusions are pro- 
vided in Section VII. 
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where the factor that renormalizes the hopping in the 
Jh —> oo limit is 



S(6i,<f>i,0j,<t>j) = cos(-^)cos(-|-) 
+ sm(-)sm(^-)e" 
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and the angles are used to parametrize the t 2 g classical 
spins at every site. The parameters , are the hopping 
amplitudes between the orbitals 7 and 7' in the direc- 
tion a. In two dimensions, these are t^ a = —V3t^ b = 

-V3t* ba = 3t£ = 1, and i*. = ^Ib = V^L = Kb = 1 
and in three dimensions one needs also consider t z aa = 
t z ab = t z ba = and tl h = 4/3. 

Although the phononic coupling is barely used in the 
effort reported here, for completeness a description if pro- 
vided. Qu, Qn and Qz% are normal modes of vibration 
that can be expressed in terms of the oxygen coordinate 
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II. MODEL AND DETAILS OF THE 
SIMULATION 



The standard double-exchange two-orbital model for 
manganites was used in this study. Since the model has 
been widely described in previous investigations^ this 
section is brief. The two orbitals arise from the e g states 
that are active at the Mn ions in Mn-oxides, as exten- 
sively discussed beforeAiiM The electronic density will 
be n=0.5 in all the results reported here. This corre- 
sponds to one electron per site of a square or cubic lattice, 
and it is the density where the CE state of manganites is 
experimentally known to be stable. 

In more detail, the Hamiltonian for the two-orbital 



Also, T xi — c l a c ib~^cl[ ) Ci a , Tzi—claCia-clfrCib, and 

Pi~c\ a Ci a +c\ b Cib- The constant A is the electron- 
phonon coupling related to the Jahn- Teller distortion 
of the MnOe octahedron^ Regarding the phononic 
stiffness, and in units of t^ a =l, the D a parameters 
arc Di~l and _D 2 =-D3=0.5, as discussed in previous 
literature^ The rest of the notation is standard, 
including a superexchange term with coupling Jaf, 
which can be made spatially inhomogeneous due to 
the introduction of quenched disorder in the form of a 
bimodal distribution of amplitude A (that will always 
be expressed in units of t^ a ), namely at each link the 
local coupling Jaf(*,o;) is either Jaf + A or Jaf — A 
with equal chance. This type of disorder is natural if the 
effect of chemical doping in manganitcs is considered, 
where ions of a given size are replaced by ions of a 
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different size creating locai stress in the crystai. Aiso, 
in some of the investigations reported here, a Zeeman 
term with field strength B was added to investigate the 
influence of magnetic fields. 

The two orbital model was here studied mainly using 
the canonical Monte Carlo method for spin-fermion sys- 
tems, widely described in previous literature £ The con- 
ductances were calculated using a method proposed by 
Verges^ and the unit of conductance throughout this pa- 
per is e 2 /h. The lattices studied here are of intermediate 
size (with the exception of those in Section VI using the 
TPEM) , but note that the computational effort was sub- 
stantial (see below) in view of the many parameters in 
the problem, plus the need of averaging over several con- 
figurations of disorder. Future work exclusively focused 
on a small region of parameter space could address even 
larger systems. 

It is also important to remark that the computational 
work described in this manuscript has several unique 
characteristics, including the vast resources that it has 
required to accurately simulate the manganite models. 
To give the reader a perspective of the effort, note that 
most of the simulations have been carried out on the 
NCCS XT3 supercomputer (with 2.4-GHz AMD Opteron 
processors and 2 GB of memory), using 100 nodes in par- 
allel for small systems and up to 1600 nodes per run on 
large lattices. The entire effort documented in this paper 
would have required several years using a cluster of com- 
modity PCs. Parallelization has been implemented with 
the MPI at three levels. First, the integration of the elec- 
tronic sector, which is carried out with the TPEM, has 
been parallelized. The next level is related to the pres- 
ence of quenched disorder in the system that requires the 
use of many different disorder configurations. Finally, the 
simulations have to sweep over a large parameter space 
(A, Jaf and temperature) allowing for a third level of 
parallelization. 



III. EFFECT OF LINK DISORDER ON THE 
RESISTANCE 

In this section, the effect of quenched disorder (in the 
Jaf link variables) on the insulating phases that appear 
in models for manganites will be investigated. Our main 
goal is to analyze if these insulators are robust, namely 
if they preserve their properties with increasing disorder, 
or whether they are fragile under the same "perturba- 
tion". Surprisingly, the latter is found to be the case: 
large changes in the properties of the insulating state oc- 
cur even for fairly modest disorder strengths. This is 
in good agreement with several experiments in the man- 
ganite context, as explained before. Our results provide 
evidence that some states, such as the CE, are fragile, 
namely easily disturbed by perturbations. The study 
in this section is organized following the sketch given in 
Fig. ^for a 3D case. First, the A-type antiferromagnctic 
state will be studied, followed by the CE phase. These 




FIG. 1: Schematic representation of the phase diagram cor- 
responding to the model Eq. Q in three dimensions^" and 
electronic density n=0.5. The influence of quenched disorder 
on the conductance will be mainly studied in this case in the 
A = plane (i. e. without electron-phonon coupling), focusing 
on the A-type and CE phases. 



are the closest states to the important FM metallic state. 
At the end of the section, the 2D CE state is also ana- 
lyzed. 



A. A-type antiferromagnet in 3D 

The A-type AF phase consists of ferromagnetically 
ordered planes, that are antiferromagnetically stacked 
along one axis. At A=0 and in 3D, this state is sta- 
ble at intermediate values of Jaf, as shown in Ref. l4pt 
It is interesting to observe that this phase is obtained 
in our study using unbiased Monte Carlo simulations^*^ 
in the clean limit (A=0) starting, for instance, from a 
fully paramagnetic (random) classical spin configuration. 
There are no problems identified here of convergences 
to higher energy states. However, note that due to its 
intrinsic anisotropy, the conductance of the A-type AF 
phase is highly directional, i. e., the system is metallic 
in the FM plane but insulating when the conductance 
is measured along the direction where the system is AF. 
In Monte Carlo simulations, while the A-type character 
of the state is always obtained in the appropriate range 
of parameters, the orientation (i. e.the antiferromagnetic 
direction) is random. This is clear since the three orien- 
tations are energetically degenerate. Since our goal is 
to study the A-type AF insulator, as it appears in ex- 
periments where it has an insulating character, then we 
measure the conductance along the insulating direction 
of the A-type phase. 

In Fig. |21 the conductance vs. temperature, T, is plot- 
ted for various disorder strengths, A. The study is car- 
ried out at Jaf = 0.15, which we found corresponds to an 
A-type antiferromagnetic ground state for A = 0. Shown 
are averages over 5 configurations of quenched disorder. 
At A=0, the conductance vanishes at low temperatures, 
result compatible with the expected insulating proper- 
ties of the A-type AF phase (the spin-spin correlations 
and spin structure factor indeed confirm the A-type spin 
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FIG. 2: Results illustrating the appearance of a metal upon 
adding quenched disorder to an insulator. Shown is the nu- 
merically calculated conductance vs. temperature T for the 
A-type AF phase in three dimensions, with the antiferromag- 
netic spin arrangement along the z— axis, with and without 
disorder. The calculation was done on a 4 3 lattice, working 
at Jaf = 0.15 and A = 0. For each temperature and A, 
5,000 thermalization and 15,000 measurement Monte Carlo 
steps through the whole lattice were carried out. The actual 
measurements of observables were separated by 10 steps. The 
starting configuration corresponded to an A-type AF in the 
^-direction. The link disorder is bimodal, Jaf ± A, and the 
average shown is over 5 quenched disorder configurations. For 
the error related to this average see Fig. QJi. 



character of this state). 

The most interesting results of Fig. [3 are obtained at 
nonzero A. Results at A=0.02 indicate an enhance- 
ment of metallicity (higher conductance) than in the 
clean A=0 limit, which is surprising. The system has 
a finite conductance in this regime. At A=0.04, (a dis- 
order strength that can be considered intermediate with 
respect to the central value of the disorder distribution, 
Jaf = 0.15), the effect reported in this paper is the most 
dramatic. At this value of A, the low temperature con- 
ductance reaches a maximum. Upon further increasing 
A, the transport properties now become that of a poor 
metal, namely the actual conductance value is low (at 
least compared with its maximum possible value), and 
the result is nearly temperature independent, ft is the 
surprising regime of small/intermediate A that it is the 
most interesting due to its dramatic change in conduc- 
tance, but the large A region is also puzzling since at very 
large disorder one would expect to find an insulator. 

In studies incorporating quenched disorder, the aver- 
ages over disorder configurations are important. Indeed, 
disorder introduces another source of error (apart from 
Monte Carlo error) and the size of this error needs to be 
determined. To that effect the dependence of the con- 
ductance on each of the five different disorder configura- 
tions used in the previous figure is shown in Fig. [3Ji for 
A = 0.1. We found an analogous behavior for other val- 
ues of A although this error tends to be larger for smaller 
As. Moreover, a similar dependence of the conductance 
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FIG. 3: (a) Dependence of the conductance on the particu- 
lar quenched-disorder configuration used (labeled as 1-5, with 
"Av." denoting the average). Although the results near the 
metal-insulator transition vary substantially from configura- 
tion to configuration, the qualitative behavior is similar in 
all of them. Shown is the conductance vs. T for the A-type 
AF phase (AF along the z-axis) for different disorder config- 
urations, using A = 0.1. The calculation was done on a 4 3 
lattice at Jaf = 0.15 and A = 0, with the same number of 
Monte Carlo steps discussed in Fig|5| The disorder is bimodal, 
Jaf ± A, and the average shown is over the 5 disorder config- 
urations, (b) Dependence of the conductance on the disorder 
strength A, at low temperatures, showing the transition from 
an insulator to a metal upon disordering the system. Shown 
is the conductance vs. A of the A-type AF phase (AF along 
the z— axis) at T = 0.02t. The lattice, couplings, Monte Carlo 
steps, initial configuration, and type of disorder used are as in 
(a). The error bars mainly arise from the average over several 
quenched disorder configurations. 



on the disorder configuration was found for the CE phase. 

The low temperature (T=0.02) dependence of the con- 
ductance on the disorder strength, A, is shown in Fig-EJa. 
It is remarkable that the connection between the insulat- 
ing state at A=0 and the large A limit is not smooth, 
but a peak is observed at intermediate disorder strength. 

A possible explanation of the results obtained in this 
study is the following. Let us start at very large A, 
namely in the regime where disorder dominates. It may 
appear unnatural that in such a limit a metal (even if 
poor) is obtained, since disorder is associated with lo- 
calization and related effects. However, in the large A 
regime each of the lattice links (defined by two nearest 
neighbor sites) has an associated pair of classical spins 
which arc cither mainly ferromagnetic or antifcrromag- 
netically aligned, depending on the value of the cou- 
pling Jaf at that particular link. If the link coupling 
is Jaf+A, then this link is AF, while for Jaf-A it tends 
to be ferromagnetic due to the double exchange mech- 
anism. The fact that at large A the orientation of the 
classical spins at each link follows the disorder is clear 
from Fig0] (top panels) where both the spin correlations 
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and the Jaf couplings arc shown. Although a few links 
have different color (or thickness) on the left and right 
top panels, the majority coincide, namely the spins fol- 
low the couplings. The important point in this regime 
is that the ferromagnetic links can percolate since there 
are (on average) 50% of them, and in three dimensions 
the critical concentration for percolation is only 25%. A 
percolative path is shown in the lower panel of Fig0] 
However, note that this explanation is only rough, since 
for Heisenberg classical variables (i. e. defined by two an- 
gles) as in our simulations the correlation between two 
links does not necessarily take two values (FM/AF) as in 
the case of Ising spins. Nevertheless, we have observed 
that even within the many percolative channels, the spins 
tend to twist their orientations to satisfy better the com- 
peting FM-AF tendencies reducing the conductance sub- 
stantially. The situation would be different in systems 
with sufficient anisotropy to be described by Ising vari- 
ables. The important fact is that even if the FM links are 
percolated, the AF links provide a frustrating tendency 
to them, i. e.they are not independent of each other. 
Thus, the resulting conductance is finite but poor. Nev- 
ertheless, the essence of the metallic behavior found at 
large A is related to the percolative picture described in 
this paragraph. 

For low and intermediate disorder strengths, the con- 
ductance is larger than at large A. A different, more 
complex explanation of the results holds in this regime, 
and this will be discussed later in this paper, after con- 
firming that the large intermediate conductance occurs 
in other systems as well. 



B. CE phase in 3D 

An investigation similar to that reported in the pre- 
vious section for the A-type AF phase was also carried 
out for the CE phase. The CE state is stabilized in the 
clean limit by simply further increasing Jaf with respect 
to the value used for the A-type AF states, as reported 
in previous investigations (see sketch Fig^l . In Fig. |SJa) 
the conductance vs. T is presented for Jaf=0.2 (corre- 
sponding to a CE-phasc ground state at A=0), and for 
various strengths of the disorder configuration. At A=0, 
the expected result is found, namely for temperatures be- 
low a critical one the conductance vanishes. This clearly 
signals an insulator, which has CE characteristics based 
on the values of the spin and charge correlations^ For 
small values of A, the conductance increases from zero 
at low temperatures, and at A=0.10 it reaches a plateau, 
which is nearly temperature independent. Note that con- 
trary to the results for the A-type AF phase, the conduc- 
tance grows monotonically until the large A (percolated) 
regime is reached. It will be argued later that this is 
caused by the fact that the FM metallic phase is no longer 
a neighbor of the CE state in the 3D phase diagram. For 
completeness, the case of a FM ground state (stabilized 
at smaller values of Jaf) is also shown in Fig.[5Ib). At 




FIG. 4: Figures illustrating the expected dominance of 
quenched disorder over the spin configurations, at large A. 
(a) Sign of Si ■ Sj for nearest-neighbor sites indicated 
as red/thin (negative sign; antiferromagnetic) and blue/thick 
(positive sign; ferromagnetic) connections. The result cor- 
responds to one particular quenched disorder and classical 
spin configurations, obtained after the Monte Carlo thermal- 
ization. (b) Value of the link coupling JAF(i,j), indicated 
as colors red/thin (Jaf + A) and blue/thick (Jaf — A). In 
both left and right panels, the couplings were Jaf=0.15 and 
A=0.1. The correlation between the results in both panels is 
apparent, indicating that for this large value of the disorder 
strength the classical spins follow the disorder in the spin-spin 
correlations, (a) An example of a FM percolative path found 
in the study of the top panels. There are many of them (but 
only one shown). 



all temperatures, the conductance decreases with increas- 
ing disorder strength. For large enough A (not shown), 
a common limit is reached both in the cases of the (a) 
and (b) panels in the figure. The metallic FM state is 
damaged by disorder and its conductance is reduced, as 
intuition suggests, while the insulating states (A or CE 
type) actually become metallic after the addition of dis- 
order. 



C. CE phase in 2D 

To verify the universality of the rapid transformation 
of the CE phase into a poor metal upon introduction 
of disorder found in 3D, similar calculations have been 
carried out in 2D. In this case, a plethora of previous 
investigations^ have shown that the CE phase is stabi- 
lized at n=0.5 by increasing the coupling Jaf- A typi- 
cal phase diagram is sketched in Fig. Note that the 
CE phase is stabilized even at A=0, mainly for topologi- 
cal reasons related with the kinetic energy of carriers in 
zigzag chains^ However, in the analysis below, we will 
also study the CE phase in the regime where the electron- 
phonon coupling A is important, to search for qualitative 
differences. 
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FIG. 5: (a) Results illustrating the metal to insulator tran- 
sition induced by disorder, focusing on the regime of the CE 
phase in a three-dimensional lattice. Shown are the conduc- 
tance vs. T of the CE-phase, with and without quenched 
disorder of strength A. The calculation was done on a 4 3 
lattice, at Jaf = 0.20, A = 0, and using 5,000 for thermaliza- 
tion and 15,000 for measurement Monte Carlo steps (the lat- 
ter actually evaluating observables at every 10 steps to avoid 
self-correlations) . The initial starting configuration had a CE 
phase with zigzag chains in the xy-plane. The disorder is 
bimodal, Jaf ± A, and shown are averages over 5 disorder 
configurations, (b) Conductance vs. T for the FM phase, 
with and without disorder. In this case, where a metal exists 
at A=0, the disorder simply decreases the conductance, as 
intuitively expected. The lattice, Monte Carlo steps, number 
of disorder configurations, and type of link disorder are as in 
(a). The superexchange coupling is Jaf=0.10. 



Figure E{a) shows the effect of disorder on the CE 
phase for various strengths A and with A = 0. The re- 
sults present interesting similarities with those obtained 
for the A- type AF phase in 3D: a conducting system is 
rapidly stabilized by increasing A and this phenomenon 
occurs in a non-monotonic way, with fairly small A's 
having the largest effect in the conductance. At the 
lowest temperatures the system has an insulating behav- 
ior that becomes metallic at intermediate temperatures 
and finally it achieves constant conductance at the high 
temperature paramagnetic regime^ It is interesting to 
observe that increasing A, the rapid transformation of 
the insulator into a metal no longer occurs, as shown in 
Fig. [7fb) for the case A = 1. While a nonzero conduc- 
tance is reached at low temperatures, its value is substan- 
tially smaller than in the A=0 case. This effect is more 
pronounced in Fig. |SJa), where the conductance at low 
temperatures T=0.01 is shown. The large difference be- 
tween the small and large A regimes is obvious: at small 
A the fragility of the CE state is much more notorious. 
Also note the similarity of the conductance vs. A plot for 
the 2D CE state with the results for the 3D A-type AF 
state. We believe the reason is that both these states are 
the closest to the FM metallic state in parameter space, 




FIG. 6: Schematic phase diagram for model Eq. Q in two 
dimensions, from RefllSl The effect of quenched disorder is 
studied in this section in the cases: (i) A = 0, Jaf=0.17, 
corresponding to a CE phase without active phonons, and 
(ii) A = 1, Jaf=0.145, corresponding to a CE phase stabilized 
including phonons. These two cases are indicated by dashed 
lines. 
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FIG. 7: Insulator to metal transition induced by quenched 
disorder in two dimensions, (a) Shown is the conductance 
vs. T of the two-dimensional CE phase, for different disorder 
strengths A. The calculation was done on an 8 x 8 lattice, 
at Jaf=0.17, and A = 0. The number of Monte Carlo steps 
was 5,000 for thermalization and 15,000 for measurements, 
calculating observables every 10 steps. The initial configu- 
ration was the perfect CE phase. The disorder is bimodal, 
Jaf ± A, and the averages shown are over 5 quenched disor- 
der configurations, (b) Same as (a) but for Jaf = 0.145 and 
A = 1.0. 



as explained in more detail below. 

Finally, the clear correlation between the spin 
orientations^ at each link and the value of Jaf at that 
particular link is shown in Fig. IH^ h.c) at large A. As in 
3D, percolative processes favor a poor metallic behavior 
in this regime. But note that 50% is the critical num- 
ber of links for 2D percolation and, as a consequence, 
the 2D large A regime is only at the verge of metallicity. 
This explains the smaller values in conductance at large 
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FIG. 8: (a) Dependence of the conductance on disorder 
strength A, at low temperature T=0.01, using results taken 
from Fig.|7Ja). The most dramatic evidence of fragility in the 
2D insulating CE state is obtained in the absence of electron- 
phonon coupling, (b) Sign of Si ■ Sj for nearest-neighbor sites 
indicated as red/thin (negative; AF) and blue/thick (pos- 
itive; FM) connections, (c) Value of Jaf(*, j) indicated as two 
colors red/thin (Jaf + A) and blue/thick (Jaf — A). In both 
cases, Jaf = 0.17, A=0.1, and the results were obtained for 
fixed classical spin and quenched disorder configurations. The 
correlation between the quantities in the (b) and (c) panels is 
apparent, indicating that for this large value of A the system 
follows the disorder in the spin-spin correlations. 



A between 2D and 3D. The state induced by the disorder 
acting on the 2D CE state must be similar to the "CE 
glass" studied in Ref. before. 



IV. EFFECT OF DISORDER ON THE 
MAGNETO-RESISTANCE 

The fragility of the insulating states in the A=0 regime 
upon the introduction of disorder, suggests that other 
perturbations may also rapidly transform these states 
from insulating to metal. In particular, the effect of mag- 
netic fields, B, is important for potential applications 
and to improve our understanding of the CMR effect of 
manganites. To carry out this study, we will proceed se- 
quentially through the phases and parameters described 
before and summarized in Figs. n an dEl The magneto- 
resistance here will be defined by [(i?(0) — R(B))/R(0)] x 
100, where R(B) is the resistance (1 /conductance) at the 



magnetic field B. The magnetic field is applied in the z 
direction, same orientation as the A and CE phases, and 
we have chosen the intensities of 5=0.05 and 5=0.10 for 
our studies. 



A. A-type AF phase 

Results for the conductance in the 3D A-type AF 
regime are given in Fig. in the clean and disordered 
cases, and considering both B=0 and 5=0.1. It is re- 
markable that in the clean limit, a relatively small field 
such as 5=0.1 is capable of totally transforming the in- 
sulator into a metal. In this case, the magneto-resistance 
factor is the largest, very close to 100%, as shown in 
Fig. ^3 A similar phenomenon occurs at 5=0.05 (not 
shown). The magnitude of this effect remains strong for 
small As such as 0.02, and then slowly decreases. The 
singular behavior at temperatures close to T=0.03-0.04 is 
caused by the peaks in conductance at those intermediate 
temperatures. The focus should be on the low tempera- 
ture regime, which is the most interesting in view of the 
magnitude of the effects. The result at A=0 is probably 
caused by the close proximity in energy of the FM metal- 
lic state, as discussed before in RefllSl If this is the case, 
then the transition must be first-order between the two, 
and a finite value of B is needed to induce it. From this 
perspective, the finite but small A region will be more 
interesting since even in the limit of vanishing 5, it has 
a substantial magnetic susceptibility, as discussed below. 

B. CE-phase in 3D and 2D 

Monte Carlo simulations with an applied magnetic field 
were also carried out in the regime of the CE phase in 
three dimensions, and the results are shown in Fig. 1111 
for B= 0.1. As in the previous case of the A-type phase, 
here also the CE phase is found to be highly suscepti- 
ble to the application of a magnetic field, and the clean 
limit insulating state quickly becomes metallic. Disor- 
dered samples are less susceptible to the application of 
a magnetic field, as previously observed for the A-type 
AF phase. Both regimes, CE and A-type AF, present a 
similar phenomenology. 

The analysis of the influence of magnetic fields can 
be made more precise by studying the magnetization at 
small fields B. The results are shown in Fig. At 
A=0, a small but finite B is needed to render the system 
ferromagnetic, as previously discussed. This is probably 
related with the first-order transition reported in previ- 
ous investigations in the clean limit>i& The slope of M vs. 
B at low fields increases with A. This small A regime 
is also the most relevant in some of the cases studied in 
previous sections in the absence of magnetic fields. 

To further investigate the origin of these effects, Monte 
Carlo thermalized spin configurations have been ana- 
lyzed. From each of these configurations we constructed 
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FIG. 9: Influence of a magnetic field on the transport prop- 
erties of the 3D A-type AF state, as discussed in the text. 
The effect is particularly dramatic at small A, including the 
clean limit A=0. Shown are the conductances vs. T of the A- 
type AF phase, at different disorder strengths A. Filled sym- 
bols (empty symbols) correspond to a magnetic field B—0.1 
(B=0.0). The calculations were done on a 4 3 lattice, at 
Jaf — 0.15 and A=0, and using 5,000 thermalization and 
15,000 measurement Monte Carlo steps (calculating the con- 
ductance every 10 steps). The starting configuration used 
was an A-type AF phase. The disorder is bimodal, Jap ± A, 
and shown are averages over 5 disorder configurations. The 
magnetic field is applied in the z-direction. 
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FIG. 11: (a) Influence of a magnetic field on the transport 
properties in the regime of the CE phase in three dimensions. 
Shown are the conductances vs. T of the CE-type phase, for 
different disorder strengths A. Filled symbols (empty sym- 
bols) correspond to magnetic fields .6=0.1 (B=0.0). Note 
the large change at low temperatures, and small and vanish- 
ing disorder strength. Lattice, A, Monte Carlo steps, number 
of quenched disorder configurations for averages, and other 
details are as in Fig|§] The coupling Jaf is 0.2. (b) Magneto- 
resistance (defined in Section HV1 corresponding to the data 
shown in (a). 
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FIG. 10: Magneto-resistance (defined in Section IIV|l at 
B=0.1 vs. T, for the A-type AF regime in 3D, and con- 
sidering several disorder strengths A (indicated). The lattice, 
couplings, Monte Carlo steps, and other conventions are as in 
Fig. 
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FIG. 12: Magnetization vs. applied magnetic field, for various 
disorder strengths A at T=0.01. Inset: The magnetic suscep- 
tibility dM / (1B\b^o vs. disorder strength A. Parameters and 
other details are as in Fig. 1111 



the "summed snapshot", X^ev ^0' where V, represents 
the vicinity (4 sites in 2D or 6 sites in 3D) of site i. The 
advantage of this quantity is that it helps identify regions 
that are FM. If the spins near i are pointing randomly 
then the "summed" spin tends to cancel. In Fig. I13f a). 
an example is shown using an 8 x 8 lattice (very simi- 
lar results were obtained on a 20 x 20 cluster using the 
TPEM method described later in this paper). In the ab- 
sence of magnetic fields, but including disorder with a 



small A, the stabilized spin configuration is very differ- 
ent from the CE phase with parallel zigzag chains that 
is present at A=0. With this disorder, ferromagnetic 
clusters containing many spins with the same orientation 
are present. Their random orientation character gives 
a nearly zero global magnetization. However, clearly a 
system with preformed FM islands must have a high sus- 
ceptibility to external magnetic fields, as already argued 
to be the reason for the manganite CMR effect in previ- 
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FIG. 13: (a) Monte Carlo "summed snapshot" showing the 
"clustered spins" stabilized by a small amount of quenched 
disorder at B=0, for the case of an 8x8 lattice, Jaf = 0.17, 
A=0, and A=0.02. In (a), the presence of FM clusters is 
clear to the eye: quenched disorder with a small value of 
A transforms the original CE state into a state with local 
ferromagnetism (but not global). The conductance of this 
state is zero, (b) Same as in (a) but for B—0.02. A very 
small magnetic field is sufficient to align the pre-formed FM 
clusters, and simultaneously increasing the conductance (to a 
value C=3, in this case). 
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FIG. 14: (a) Conductance vs. T of the two-dimensional CE 
phase for different disorder strengths A. Full (empty) symbols 
correspond to a magnetic field B=0.1 (B—0). The calculation 
was done on an 8 x 8 lattice at Jaf = 0.17, A=0, and using 
5,000 thermalization and 15,000 measurement Monte Carlo 
steps, and the CE phase as initial configuration . The disorder 
is bimodal, Jaf±A, and the average shown is over 5 quenched 
disorder configurations, (b) Magneto-resistance at B=0.1 vs. 
T for the 2D CE phase, at different disorder strengths A, as 
indicated. The calculation was done using the same lattice, 
couplings, and other details as in (a). 



cms studies^ When magnetic fields are applied, the FM 
clusters rapidly orient their moments creating a global 
FM state, as shown in Fig. I13f b1. 

For completeness, the effect of a magnetic field on the 
CE phase in two dimensions was also studied. Results for 
the conductance and the magnetoresistance are in Fig. 1141 
for A=0. The trends arc very similar to those found in 
three dimensions. 
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FIG. 15: (a) Resistance vs. temperature for a 12x12 lattice 
in the clean limit A=0, A = 0, and considering Jaf = 0.15 
(squares), Jaf = 0.16 (circles), and Jaf = 0.17 (triangles). 
The figure illustrates the metallic character of the Jaf = 0.15 
state (with a FM ground state), and the insulating proper- 
ties of the Jaf = 0.16 and 0.17 samples (with a CE ground 
state). The Monte Carlo simulation was carried out starting 
with a random state at each temperature, using a warming 
up of 10,000 steps, followed by an additional 10,000 steps 
to measure properties, (b) Same as (a) at Jaf = 0.17, but 
now considering the on-site disorder strengths Ap indicated, 
and averaging over 3 different disorder configurations. The 
resistance of the clean sample is reduced several orders of 
magnitude by the addition of the on-site potential disorder, 
similarly as it occurs for link disorder, (c) Same as (a), but 
on a 4 3 lattice and with Jaf = 0.20. 



V. RESULTS USING ON-SITE DISORDER 

In this section the influence of an on-site disorder (bi- 
modal, of strength Ap) will be investigated. The term 
added to the Hamiltonian is an on-site random energy, 
Y^i 7 £ic| 7 Ci ;7 with £j = ±Ap, to be contrasted with the 
link randomness considered in the previous sections. It 
will be shown that the conclusions regarding the fragility 
of the insulating states at A=0 remain valid, even if the 
source of disorder is changed. 

Figure fTST a) shows the resistance obtained in the clean 
limit, parametric with Jaf, using a 2D cluster. This re- 
sult is the analog of the conductance plots in the clean 
limit presented in previous sections. Varying Jaf, &n 
abrupt transition from a metal to an insulator is ob- 
tained at low temperatures, as expected. Let us now 
introduce on-site disorder. The results are in Fig. I15f b). 
It is clear that the insulating state dramatically reduces 
its resistance with increasing Ap , transforming into a bad 
metallic state, as found in the case of link disorder. The 
changes in the resistance are of several orders of magni- 
tude. The same occurs in three dimensions, Fig. H5lV ). 
where in addition it is shown that the resistance has a 
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FIG. 16: (a) Effect of quenched disorder on the two- 
dimensional CE phase, running the simulation with the 
TPEM. Conductance vs. T of the CE phase on a 20x20 
lattice calculated with the TPEM method using M=30 mo- 
ments, and the TPEM parameters e pr =10~ 6 and etr=10 -6 . 
The coupling was Jaf = 0.17, density n=0.5, and 2,000 for 
thermalization and 1,000 for measurement Monte Carlo steps 
were carried out. The CE phase shows zero conductance at 
low temperature in the clean limit. This transforms into a 
disordered metal with finite conductance when link quenched 
disorder of strength A=0.05 is added, (b) Effect of disor- 
der on the three-dimensional CE phase, studied with TPEM. 
Details are the same as in (a) but on an 8 3 lattice and at 
Jaf = 0.20. 



minimum at intermediate values of the disorder strength 
Ap, as it occurs for link disorder. Although the substan- 
tial numerical effort involved in the simulations presented 
in this manuscript prevents us from discussing the site 
disorder case in as much detailed as the case of link dis- 
order, nevertheless it is safe to conclude that site disorder 
has a phenomenology similar to that of link disorder, in 
the sense that by increasing Ap a conductor is obtained. 



VI. FINITE-SIZE EFFECTS AND TPEM 
TECHNIQUE 

The investigations reported in the previous sections 
have been carried out using intermediate size clusters, at 
the verge of the sizes that can be handled with the exact 
diagonalization treatment of the fcrmionic sector. Fortu- 
nately, the systematic effects observed, and similarities 
of results between different quenched disorder configu- 
rations, dimensionality, and insulating phases strongly 
point toward a robust effect, and we are confident that 
the essence of the phenomenon has been unveiled. Nev- 
ertheless, only the study of larger clusters would con- 
vincingly rule out spurious size effects in our investiga- 
tions. Fortunately, recent algorithm developments have 
led to the use of the "Truncated Polynomial Expansion 
Method" (TPEM) 43 i 44 i 45 i 46 as an alternative to conven- 



tional Monte Carlo methods where the fermionic sector 
is treated exactly. The advantage of the TPEM is that 
its CPU time grows linearly with the number of sites N, 
while the canonical approach used in the previous sec- 
tions demands an effort that scales as A^ 4 . This advan- 
tage is not without penalization: the coefficient in front 
of N in the TPEM scaling is very large. For this reason, 
we studied few large lattices, mainly to check if the re- 
sults on small systems are preserved at larger N values. 
After confirming that the TPEM reproduces the results 
of small systems accurately (not shown but see Ref . , 
studies were carried out using 20 x 20 and 8x8x8 clus- 
ters. In Fig. the results are shown, both in the clean 
limit and introducing disorder. In both cases, at low tem- 
peratures the conductance changes from exactly zero to 
a finite number, as observed in previous sections. While 
this is not a final proof, at least it has been confirmed 
that increasing the lattice sizes does not lead to dramatic 
changes in our conclusions, and the insulator to metal 
transition with increasing A seems to be a solid property 
of the CE states of manganites at small electron-phonon 
coupling. 



VII. DISCUSSION AND CONCLUSIONS 

In this paper, the influence of small amounts of 
quenched disorder on the transport properties of man- 
ganite models was investigated. The focus has been on 
the insulating phases that are stabilized without electron- 
phonon coupling)** namely those that have insulating 
properties due to geometrical considerations and frustrat- 
ing competing tendencies between FM and AF states. 

The study reported here has identified regimes where 
the insulating phases, the CE and A-type AF states, are 
very susceptible to the introduction of disorder. For in- 
stance, the A-type AF phase in 3D and the CE-type 
phase in 2D, they both present a conductance that 
sharply increases with increasing A from zero, and then 
it settles to a smaller value at large A. The reason ap- 
pears to be explained with the type of spin arrangements 
that become stabilized in the presence of quenched disor- 
der. Typical results are shown in Fig. El f° r the case of 
the 2D CE state. In (a), the CE pattern is nicely repro- 
duced by the Monte Carlo simulation at zero and very 
small A. However, in (b), a key observation is that lines 
with classical spins align ferromagnetically are stabilized 
along the vertical direction for a small A. For instance, 
the first and fourth column from the left are all FM. 
The second, fourth and fifth columns are almost entirely 
FM, with the exception of only one link. More quanti- 
tatively, there are 50 links FM in the vertical direction 
but only 27 FM links in the horizontal direction (ratio 
approximately 2). These lines provide conducting paths 
for the charge, rendering the system conducting in the 
direction of those lines. Even though the conductance is 
not always measured along these lines, the average over 
all directions is still large. This is to be contrasted with 



12 



J:,\ 









































































































































Hh 

















































































































FIG. 17: Results of Monte Carlo simulations that summarize 
the main conclusions of this manuscript (see text for discus- 
sion). Shown are the signs of the classical spin correlation 
Si ■ Sj considering nearest-neighbor sites i,j, indicated as 
red/thin (negative; AF) and blue/thick (positive; FM) con- 
nections. The lattice is 8x8, Jaf = 0.17, and the disorder 
strength is (a) A=0.01, (b) A=0.02, (c) A=0.08. In (a), the 
disorder is too weak to change the CE spin pattern, which is 
clear in the result. However, in (b) the CE phase is replaced 
by long ferromagnetic (metallic) chains, which in this case 
are vertical. These chains provide conductance channels in 
the vertical direction. An equal number of horizontally con- 
ducting systems are found when analyzing many quenched 
disorder realizations, thus restoring rotational invariance af- 
ter averaging. In (c), the state is more reminiscent of the high 
disorder regime, with shorter chains. At large A, the spins 
follow the disorder as explained already in Section fill CI Of 
the three configurations, the largest conductance is obtained 
in (b). 



the regime shown in (c) for large values of A, where the 
ratio between FM links along the vertical and horizontal 
axes is approximately 1, as is also in the clean limit. At 
least for 3D, this type of configurations is also conduct- 
ing, as explained with a percolative picture in previous 
sections. These tendencies were observed in numerous 
realizations of disorder, with the only difference between 
them being that the dominance of the horizontal or ver- 
tical lines has equal chance. Although most of the inves- 
tigations reported here were obtained using intermediate 
size clusters, some larger systems were studied with the 
TPEM and the results remain similar. Nevertheless, a 
more sophisticated analysis is, of course, always desir- 
able to confirm our conjectures. Hopefully, our effort will 
induce further work in this challenging area of research. 

The origin of the FM lines that arc rapidly stabilized 
by disorder must be the close proximity in energy of the 
full FM state. While the CE state prefers zigzag chains 
(which are band insulators^) to minimize the energy, 
the FM metallic state prefers straight lines since they 
are metallic. So this striped state can be considered as 
a compromise between the CE zigzag insulating and the 
FM uniform metallic states of the clean limit. Note that 
it is also possible to imagine a "clustered" state as a 
compromise of these two tendencies, and indeed in some 
circumstances this type of states was observed in our sim- 
ulations, as it was shown in Fig. ll3f a). But this state has 



zero conductance, while the striped state of Fig. lTTT b) has 
a finite conductance and, thus, it is responsible for the 
rapid insulator to metal transition with increasing A. 

We believe the conclusions reached in this paper have 
consequences beyond the case of Mn oxides. Clearly, 
the entire physics discussed arises from competition 
of phases, which is a common phenomenon in transi- 
tion metal oxides. This phase competition is already 
present in the clean limit, since the CE and A- type AF 
phases appear in parameter space in narrow regions of 
the phase diagram, as a compromise between the FM 
and AF dominant tendencies. The same was recently 
shown to occur in between a d-wave superconductor and 
an antiferromagnetSi in a study of a phenomenological 
model for cuprates: in this case a striped phase was found 
as a compromise of the two tendencies in some regions 
of parameter space. It should not be surprising that the 
intermediate phases stabilized by phase competition arc 
fragile, namely they change dramatically when influenced 
by disorder or external fields. This fragility manifests 
in the fairly exotic states that are stabilized when small 
quenched disorder is introduced: for the case of the man- 
ganites, these are states with either FM clusters or FM 
lines as reported here, having in common regions where 
ferromagnetism and mctallicity dominates locally. That 
this occurs at values of A as small as 0.02 is the most sur- 
prising result. This high susceptibility to disorder should 
exist in many other complex oxides as well. When these 
giant effects are compared with other insulator to metal 
transitions induced by disorder (see Section it appears 
that the cases discussed here have a unique character, not 
previously discussed in the literature. The consequences 
of these results for our understanding of manganites and 
the CMR effect seem clear: it is the fragile nature of the 
insulator to the introduction of external perturbations 
(quenched disorder, magnetic fields, electric fields, . . . ) 
that may cause the well-known giant responses of Mn 
oxides. 
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